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Abstract 

We address the strong-coupling regime effects in the light-like quark jets where the ra- 
diated gluons are hard but highly collinear. These may lead to additional contributions 
in the invariant mass (or recoil energy) spectra, on top of the Fermi motion of the initial 
heavy quark in decays like b— > s + 7. It is shown that the integer moments are nevertheless 
free from such effects; perturbation theory corrections for the moments are driven by a s 
in the small coupling regime. They are modified nonperturbatively by the soft modes but 
not by the collinear modes. 
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1 Introduction 



The OPE provides a consistent approach to describe sufficiently inclusive distributions in 
the decays of heavy flavor hadrons. It incorporates power-suppressed effects which orig- 
inate from the nonperturbative domain that is responsible for confinement and bound- 
state dynamics. Its important ingredient is the separation between the contributions of 
different momentum scales, with the short distances generically referred to the coefficient 
functions, and large distances ~ Aqq D belonging to the nonperturbative operators. This 
general approach applies to both the weak transitions into sufficiently heavy quarks like 
b^ciu, and to heavy-to-light decays like b^s + '-f or b^uiv. The former case has been 
best elaborated including the technical implementation of the Wilsonian 'hard' momen- 
tum scale separation in the perturbative expansion. This approach appeared successful in 
describing experimental data [1]. There remain, however certain principal differences be- 
tween the decays into heavy and light-like quarks when gluon bremsstrahlung is carefully 
incorporated. 

The standard application of the OPE is initially formulated in the Euclidean theory. 
It utilizes the expansion of the amplitudes in the momentum of the soft quantum fields 
which describe nonperturbative dynamics. Therefore, it assumes all the components of 
the momentum to be bounded by some mass \x identified with the normalization scale of 
the nonperturbative operators. It is normally chosen to be 



where /ihadr — 600 MeV represents the typical momentum scale characteristic for the strong 
coupling domain in QCD. Consequently, for the gluon exchanges those with any compo- 
nent of the momentum k a larger than [i are attributed to the coefficient functions and are 
calculated in practice in perturbation theory. The same idea is transferred with minimal 
modifications to decays in the Minkowski kinematics. 

This procedure works fine for the transitions between heavy quarks, in particular in 
b — > civ. The lower cutoff on the gluon energy automatically eliminates the infrared 
domain and no running coupling a s (q 2 ) enters at the scales q 2 <C /i 2 . The perturbative 
expansion is therefore applicable and justified. 

The situation looks different in the heavy-to-light decays, in particular in b — > s + 7. 
The gluon bremsstrahlung is described by the double logarithmic probability 



where dWbom is the bare 'hard' decay rate and = (u, kj_, k\\ ) is the gluon four-momentum. 

fen 

This expression literally assumes cu <C m&, k± <C and the radiation angle 9 ~ -J- <C 1. 
More accurate expressions for arbitrary 9 do not change the principal point here. 

The radiation probability (1) underlies the conceptual problem: even if the gluon is 
very energetic by itself, 00 ^> /i, the transverse momentum k± can be low, \k±\ ^ Aqcd 
if the gluon is highly collinear with 9 <C ^ < 1. The 'perturbative' bremsstrahlung then 



fi ~ few x /i hadr , 
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runs into the strong-coupling domain and cannot be unambiguously evaluated expanding 
in a s (Q 2 ) with Q 2 > /i 2 . This observation may cast doubts on calculability of even the 
'robust' observables like the moments of the decay distributions. 

A more careful analysis, however suggests that the primary OPE results for the (in- 
teger) moments relating them to the local heavy quark operators, remain unchanged, 
while the spectra themselves, in general, may possibly depend to some extent on new 
strong-coupling effects appearing in jet physics. 



2 OPE for the decay distributions in b^q 

For the sake of simplicity we will phrase the subsequent discussion for the decay b — > s + 7 
(or, generically, Q — > q + ip with m q = and m v = 0, referring to the colorless ip as a 
photon whether or not it carries spin). We will also assume a much stronger hierarchy 
fnb = m Q /^hadr than exists in reality, so that we may discard inessential l/nib corrections 
ab initio. 

The photon spectrum in b — > s + 7 is represented by the convolution of the soft 
'primordial' distribution function F(k + ) often known [2, 3] as describing "Fermi motion" 
[4, 5] of the b quark inside B meson, and of the 'hard' spectrum evolving from the bare 
two-body S(E 1 — ^y-) due to the gluon bremsstrahlung: 

= j_Jk + F(k + ) ^ (£,-*) • (2) 

The support of F(k + ) actually lies below A ~ Mb — rrib(n); at large negative arguments 
with I k + 1 ^> yUhadr the distribution function F(k + ) must decrease exponentially. 
Kinematically we have 

Ml -Ml 

R = x , 3 

7 2M B 

where M\ is the invariant hadronic mass squared in the final state (jet invariant mass, in 
the perturbative description). Instead of the photon spectrum we can, therefore speak of 
the distribution in Mj-, or of the moments of the hadronic mass squared which has more 
universal application in jet physics. The typical nonperturbative (bound-state) domain 
in M\ (referred to as 'window' in Ref. [3]) is M|-^Am fe , and larger M\ emerges due to 
hard bremsstrahlung. 

As schematically illustrated in Fig. 1, F(k + ) describes the intrinsic properties of the 
decaying bound state. It is universal with respect to any concrete type of Q^q transition 
as long as the recoiling system is colorless and the energy of the light quark is large. All 
field modes included in it have wavelengths limited by /i. Therefore, it is completely 
independent of mj once power corrections ~/ihadr/^b are neglected. 

F(k + ) is constrained by its (integer) moments like the usual DIS structure functions. 
The moments are given by the expectation values of the local heavy quark operators over 
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Figure 1: The corrections to the inclusive decay distribution in the OPE illustrated with the 
perturbative diagrams. The red gluon line shows the 'soft' modes with all the components of 
the four-momentum small. They are responsible for the 'Fermi motion' and are associated with 
the decaying bound state. The blue gluon line stands for harder gluons with u > [i attributed 
to r per t, which nevertheless can be highly collinear and emitted with large coupling. 



the initial hadron: 
^o = 
F 2 = 
F n = 



dk. 
dk. 
dk. 



F{k + ) - 
k\ F(k A 
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2M, 



Fi = J dk + k + F(k+) = , 
F 3 = Jdk + k 3 + F(k + ) = -^, 
iH Q \QiD z {iD,-iD z ) n - 2 iD z Q\H Q ) . 



(4) 



In Eq. (2) 



dr pcrt 



accounts for all other gluon modes. In practical terms, the renor- 



malization scheme [6] is advantageous where, at the level of one gluon emissions the 
normalization procedure sets the separation based on the gluon energy u. Then, for ra- 
diating a single gluon we have to perform the integration with only the lower cutoff on u, 
for instance 



dM 2 x 



duo 



dk 2 



Cj 



a s (k] 



Ti 



S(M 2 x -k 2 ± ^ 



(5) 



Here m&/2 simply stands for the initial energy of the s quark jet Ej Ct . The kinematic 
enhancement factor m&/(2o;) in M\ brings in the contribution of large transverse distances 
x± > /ihadr ^ ne calculation of the perturbative spectrum in Eq. (5) once M\ becomes 
smaller than ^ /ihadr^^ 1 • The factor is smaller than unity, yet it represents a numeric 
rather than parametric suppression. Perturbation theory alone is unable to precisely 
calculate its decay spectrum closer to the endpoint than a fraction of the distribution 
function domain (window), /ihadr^V This fraction appears to depend on the normalization 
point. Consequently, the full point-to-point spectrum may not be completely expressed 
through the distribution function with the resolution arbitrarily higher than the size of 
the Fermi motion 'window'. 
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It is important that the type of the strong-coupling effects generated by the 'pertur- 
bative' jet function Eq. (5) is physically distinct from the Fermi motion. The latter is 
determined by the bound-state dynamics and depends on the B meson wavefunction at 
distances ~ l^l&r from the heavy quark and its decay vertex. The jet hadronization, on 
the contrary, runs into the domain of the nonperturbative coupling only over the time 
interval ~ A*hadr m the f rame accompanying the final light quark, a much longer period 
r ict ^ /^haHrA/^due ^° the Lorentz slowdown. The hadronizing system travels over a 

1 _ 3 

distance > m fe 2 u ha ^ r before that. The nonperturbative hadronization occurs far away in 
space from the decay point where the original B meson was located. This suggest that 
the considered nonperturbative jet effects are largely universal and are independent of 
the properties of the initial bound state. There are arguments that this separation is 
exponential in the ratio /V/ihadr- 1 We will return to this point later. 



3 Moments of the distribution 

The incalculable nonperturbative end point in and the associated indeterminacy of 

%tjt raise the concern of whether the moments of the total spectrum can still be obtained, 
in the presence of gluon bremsstrahlung, in terms of the local heavy quark expectation 
values applying only the truly short-distance corrections evaluated in perturbation the- 
ory over the small-coupling domain. At first glance, Eq. (5) mandates the presence of 
additional nonperturbative corrections from jet physics at small kj_, at least at the level 
(mb/ihadr) n -(/-thadr//-t) ri for the nth moment. In the zeroth moment of , the total decay 
rate r per t such effects cancel between the bremsstrahlung and the virtual corrections as 
required by the KLN theorem. The virtual effects are absent from the positive moments, 
and the simple-minded application of Eq. (5) yields non-vanishing nonperturbative con- 
tributions already to (M^) proportional to the first moment A2 = —A 2 of the effective 
coupling (the log-moment of the dispersion coupling a e s(s)) [7] parameterizing power 
corrections in a number of observables in jet physics. 

We find, however that such a conclusion would be incorrect. The nonperturbative jet 
effects, although possibly present in the differential distribution itself, cancel out in the 
integer moments. The bremsstrahlung corrections to the usual nonperturbative expres- 
sions for the moments in terms of the £?-meson local heavy quark expectation values, 
calculated in perturbation theory, do not involve the running coupling at scales below u. 2 

The most transparent way to illustrate physics behind the related 'conspiracy' in the 
strong-coupling jet effects uses the so-called dispersive approach to perturbation theory 
in Minkowski space scrutinized by Dokshitzer, Marchesini and Webber [8]; a recent brief 
summary for purely perturbative BLM-type applications can be found in Refs. [9]. The 

: Yu. Dokshitzer, private communication. 

2 We do not consider the moment rank n as a large parameter and, therefore do not distinguish between, 
say a s {n) and a s {^/^/n). 
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method is based on the dispersion representation for the dressed gluon propagator 



af(Q 2 ) fdX 2 p(A 2 ) 1 

via a fictitious gluon mass A; Q 2 is the Euclidean virtuality and af denotes the stan- 
dard Euclidean running coupling (the above expression for simplicity assumes no pole in 
af(Q 2 )/Q 2 at Q 2 = 0). This requires calculating the inclusive observable in question O to 
order a s with an arbitrary gluon mass A, O hoin + C F ^ Oi(A 2 ). The effect of the running 
a s is then obtained by 

resum = Q hom + Cp j _ p ^ ^ 

(the order- (3 a 2 s version of this BLM improvement was formulated in Ref. [10]). In this 
context the OPE separates the effects into 'perturbative' and 'nonperturbative' according 
to their behavior in the deep Euclidean domain, not for the Minkowski objects. Therefore, 
one puts 

p(A 2 ) = p pcrt (A 2 ) + MA 2 ) (8) 

and it is assumed for the corresponding dispersion integrals, af (Q 2 ) and 5af(Q 2 ), that 
the latter dies out fast at large Q 2 > 0. 

Being interested in the effects from the strong-coupling regime we need to retain only 
the Saf (Q 2 ), or 5p(\ 2 ) piece. For simplicity we will use the full coupling in the following 
equations, however, implying that the 'nonperturbative' part can be separately considered 
when required. 

The standard bremsstrahlung probability within this framework takes the following 
form: 

9 a (k 2 ) 

The integral over A is nothing but the dispersion representation for . 

K ± 

^ p(A) fciT^-^^F' (10) 

justifying the standard prescription Eq. (1). 

As already mentioned, there is the KLN cancellation of soft gluons between the 
bremsstrahlung and the virtual corrections. The latter do not contribute to the higher 

moments, however we get generally the power-divergent integrals like J^^(k 2 L ) n a s (k 2 L ). 

The effective a s (A; 2 L ) then has to be accounted for with higher accuracy. To do this we 
write explicitly, e.g. for the first moment 

(Mi)P- = C F J AM 2 X /tr^) J^% M x 6 ( M x-( k l+ X2 m) ■ (11) 
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Performing the integration over M\ yields 



The integral has split into the product of independent integrals over k\ and over A 2 , 
which means that the effective coupling enters only at the high scale ~ y/jlrfib or higher, 
cf. Eqs. (14), (16), determined by the effective ultraviolet cutoff in the integral. Formally, 
the absence of the nonperturbative contributions follows from 

r A \2 

* i "P" 6P(X2) = Q&J afiQ2) = • (13) 

The expression for the spectrum itself helps to interpret this general result: 

dr pcrt duJ dA 2 dk 2 



Ml 



This differs from the expression 



M 



2 



^L^-fj,) sy ™» x> (15) 

J CO IX 



which would literally follow from the prescription (5). The effective coupling in the 
spectrum <5 S (^M|-) differs from the usual Euclidean coupling a s (^M^): 

rQ' 2 c\\ 2 r°° HA 2 

a s (Q 2 ) = ^J o ^P(A 2 ) while af(Q 2 ) = 7rj o p(A 2 ) (16) 

The two couplings agree in the perturbative domain with 'log' accuracy where the 
momenta appear ordered in their scale. Yet they become different in the strong coupling 
regime. In particular, the moments of 5af(Q 2 ) are nonzero and, for the simplest ansatse 
[8] all have the same sign scaling as /i^dr- The integer moments of Sa s (Q 2 ), on the 
contrary, all vanish [11, 8] according to the generalization of the relation (13): 

fr\\ 2 

n j — \ 2n 5p(\ 2 ) = (-ir Q \im oo Q 2n 6afm- (17) 

While the general expression for the probability of the gluon bremsstrahlung, Eq. (1) 
remains valid, in a sense, even nonperturbatively, the actual physical effect differs from 
what would follow from Eq. (1) at the power level. The latter therefore has limited 
applicability. The reason is that at this level the one-particle massless-gluon description of 
the interaction becomes incompatible with running of a s , a field-theory effect associated 
with a few particle states of somewhat different kinematics. Having the gluon split, 



6 



while not modifying the total rate, changes fine kinematic details thus reshuffling the 
distributions. This modifies the effective coupling for differential distributions, making it 
observable-dependent once the result is forcibly interpreted in terms of the massless one- 
gluon framework. Additionally, these effective couplings obey certain integral constraints 
which effectively recognize that the Euclidean image of the inclusive bremsstrahlung is a 
genuinely short- distance process. 

The result of the calculations supports the interpretation that, in the hard collinear jet 
configurations the growth of a s from the initial a s (E^ et ) to a s (k±) is an effect of the final- 
state interaction, viz. jet splitting. As such, it is not expected to affect fully inclusive truly 
short- distance characteristics like the total decay rate or recoil moments. We see that the 
long-time jet contributions disappear in just the integer moments of M|, indicating that 
M\ is the right kinematic variable. 

4 Discussions 

Jet hadronization effects may affect the hadronic mass distribution and, consequently, the 
point-to-point recoil spectrum in the decays like .B— >X, + 7 or B^X u £u. Yet it is shown 
that there remain no such nonperturbative strong-coupling-domain contributions to the 
moments. The perturbative corrections to the OPE relations for the moments come from 
the small coupling regime. Nonperturbatively the latter are shaped by the soft modes but 
not by the collinear modes. 

There appears an interesting analogy between the strong-coupling effects in high- 
energy jets and local duality violation in the heavy quark decays specifically studied in 
Ref. [12] in the instanton vacuum ansatz. In both cases the nonperturbative effects were 
present and enhanced close to the end point in the differential distributions. Integrating 
them over the available kinematic domain might seem to only increase the effect. Yet for 
the right moments the integrated effect is always decreased being driven by the suppression 
in the corner of the underintegrated domain; it is minimal for totally inclusive rates. This 
similarity, perhaps has its roots in the general properties of the OPE, although the precise 
structure of the latter for jet physics still remains to be understood. 

The identified property of the large-energy low-k± gluon radiation leading to the ab- 
sence of the soft domain contribution from the integer moments, seems to have something 
in common with the observation by Beneke and Braun [13] of the vanishing of a certain 
class of the leading, linear in Aqcb/E corrections in the Drell-Yan production. 

The literal incalculability of the spectrum itself in the large-mQ limit near the end 
point may still seem unsatisfactory. It may be that there exists a further conspiracy in 
the small-Zc^ jet splitting which eliminates this. A possible source for insights is to examine 
the formal /z-dependence of the incalculable piece. The spectrum must be \i- independent; 
at the same time, the kinematic M^-domain affected by large-x^ physics appears to 
depend on the lower gluon energy cutoff /i. The properties of the thus defined F(k + ; fi) 
have not been studied in detail, however, which precludes us from definite conclusions. 
The problem deserves further dedicated analysis. 
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It is conceivable that such nontrivial effects at the interface of the gluon bremsstrahlung 
and of the genuine nonperturbative corrections in the point-to-point spectrum are related 
to the specific way one defines the primordial distribution function assuming the Wilsonian 
cutoff on all the components of the momentum of the soft modes. Such a definition is 
motivated by the physics of the heavy quark bound state, but it looks foreign to the 
light cone approach. In the infinite momentum frame accompanying the jet the light-cone 
k± <-> x± and E + ~ M\ are natural variables. Using them for the separation of mass 
scales would allow, however large wave vectors k\\ which are foreign to the bound-state 
dynamics. 

Including the \ow-k± high-energy jet hadronization effects into the definition of an 
effective distribution function F{k + ) may seem to be a practical way to get around the 
incalculability of the perturbative spectrum. This idea might seem to be supported by 
the mentioned universality of the jet hadronization. The advantage is doubtful, however. 
Long-wavelength 'Fermi motion' and highly collinear jet interactions seem to describe 
quite distinct physics, as mentioned in Sect. 2. Likewise the jet effects have little in com- 
mon with the B meson expectation values of the local heavy quark operators determining 
the moments of F(k + ), or with the heavy quark bound-state dynamics in general. More- 
over, they would a priori bring in the intrinsic (logarithmic) dependence of F{k + ) on the 
high scale uiq, from which the conventional definition is free. 

The strong-coupling regime effects addressed above are not specific to the weak decays 
of heavy flavors, but are common to general high-energy light-like jet processes in hadronic 
physics. We specifically phrased the discussion in the context of the inclusive heavy quark 
decays since here it can be put in a more rigorous context of the OPE, absorbing the soft 
modes into the heavy quark distribution function F(k + ; ji). 

In our analysis no indications were found towards the specific additional intermediate 
momentum scale ~Aqcd^& for the inclusive decay distributions. The bremsstrahlung in- 
tegration runs over all virtualities: the concrete domain contributing is simply determined 
by the energy of the radiated gluon, ranging from small virtualities up to m^. 
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